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SUMMARY 


In steering an automobile the driver must basically control the direction of the 
car's trajectory (heading angle) and the lateral deviation of the car relative to a 
delineated pathway. This paper considers a previously published linear control model 
of driver steering behavior which is analyzed from a stability point of view. A 
simple app’-oximate expression for a stability parameter, phase margin, is derived in 
terms of various driver and vehicle control parameters, and boundaries for stability 
are discussed. 

A field test study is reviewed that includes the measurement of driver steering 
control parameters. Phase margins derived for a range of vehicle characteristics are 
found to be generally consistent with known adaptive properties of the human opera- 
tor. The implications of these results are discussed in terms of driver adaptive 
behavior . 


INTRODUCTION 


Analysis of the closed- loop dynamic behavior of the driver/vehicle system can give 
some insight into driver and model behavior, and provide simplified analytical expres- 
sions for relationships between various model parameters. Here we will" use a linear, 
two degree of freedom vehicle model. The two degrees of freedom to be considered are 
heading, or dir#>ction of vehicle motion, and lateral position. These two variables 
are under direct control of the driver. A third basic vehicle mode not considered 
here is roll angle. Although roll angle may influence driver behavior, it is nor con- 
trolled per se by the driver. The discussion here will include a summary of previous- 
ly published work (References 1 and 2). 

The two degree of freedom model should be considered as an equivalent or approxi- 
mation to higher degree models. Thus it subsumes such things as roll steer and weight 
tr.msfer effects to a first approximation, so that the heading response of the model 
is an adequate approximation to a real vehicle for similar inputs. 


DRIVER/VEHICLE SYSTEM MODEL 


A block diagram of the driver/vehicle system model is shown in Figure 1 . The 
vehicle equations generate side velocity (v) and yaw (heading) rate as a function of 
steering inputs through the and G\ transfer functions, respectively. Kinematic 
equations then compute vehicle heading angle and lateral lane position from side 
velocity and yaw rate inputs. The driver finally develops steering corrections based 
on perceived heading and lane position errors as processed by the behavioral transfer 
functions and For the closed-loop analysis reviewed here we will consider 
steering against disturbances applied at the steering point as shown in Figure 1 . 


*This work was partially funded by the Automotive Safety Affairs Office of the 
Ford Motor Company. However, the contents of this paper represent the views of the 
author and do not necessarily reflect the official views or policy of the Ford Motor 
Company . 


‘At this point a question might be raised as to why the Y. block is not placed in 
the pathway but rather is in the pathway, which is the sum of the heading error 
and some function Vy of lane position" error y^. This arrangement is consistent with 
the perceptual info’^rmation most readily available to the driver, which is further 
described in Ref. 3* 
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Figure 1 . Driver/Vehicle Model with Two Degree of Freedom Vehicle Dynamics 


This disturbance input will be used to approximate the lateral response effects of 
roadway disturbances on the wheels, and various forces and moments caused by roadway 
slope and lateral wind gusts. 

The dynamics and stability of the Figure 1 system can best be analyzed by con- 
sidering a steering disturbance signal (6^ as the input to a closed-loop system, and 
then analyzing the total open- loop transfer function between 6^ passing through the 
vehicle dynamics and driver behavior to the point. In this way we subsume the 
multi-path portions of the system as we progress from the single control input varia- 
ble to tne single control output variaole 6^. 

Through simple block diagram algebra we can now derive an expression for the open- 
loop transfer function 6^/5^ as follows. First the heading angle and lane position 
errors cn be expressed in response to 6^ inputs as 

♦e - ; >'e ■ r > 

Next steering response can be expressed in terms of perceived heading angle and 
lane position errors as 


^yye) (2) 


Combiing Equations 1 and 2 we can then express the total open-loop transfer function 
as 



This expression can be further simplified if we now express the transfer function 
between heading and lane position to control input in terms of the vehicle dynamics 
and kinematic equations: 
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Combining Equations 3 and 4 we have 
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Now consider models for each of the component transfer functions in Equation 5. 


Vehicle Dynamics 


Two degree of freedom vehicle dynamics have previously been analyzed in some 
detail (Reference 1) from which the following material has been summarized. In gen- 
eral Che transfer functions between heading angle and lane position to steering con- 
trol input can be expressed by second-order equations. 
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where T- is the basic time constant of the vehicle's heading response and the remain- 
ing coefficients are functions of various vehicle parameters. 

Using Che few basic vehicle parameters described in Figure 2 and some simple 
assumptions we can now express Equations 6 and 7 in terms of vehicle characteristics 
as follows. The inverse of the heading time constant is given by 
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If we now assume that Che vehicle radius of gyration (k.) is approximately equal to 
the geometric mean of the axle to c.g. distance: 

•'z ■ —■ - (9) 



Cocff (Yg^) Yaw lr.ert»a{Ij2) Coeff(Yg,) 
Radius of 
Gyrotion(I^^/m) 


Figure 2. Vehicle Parameters for Steering Dynamics 
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and incrcdace two additional parameters, the stability factor K (sec / ft^), which is 
related to the SAE understeer/oversteer gradient (Reference 1) 


K (deg/sec) 


and the "axle load ratio" 


V 


we can then write approximate expressions for the remaining coefficients in Equations 
6 and 7. 
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At best, the above expressions are good approximations for many cars. At worst 
the expressions should give us a qualitative feel for the lateral dynamic characteris- 
tics that are of importance to the driver. 

Inspection of the above equations can give us insight into vehicle dynamic re- 
sponse characteristics that are important from a driver control point of view. It is 
obvious that the heading time constant (T^) dominates the vehicle dynamics, and that 
it is a direct function of speed (Equation 8). In Equations 6 and 7 the numerator and 
denominator roots are an inverse function of speea (i.e., ) , so as the vehicle 
increases in speed the heading response becomes slower (lower frequency). The sta- 
bility factor K can also exert further influence as a function of speed. For an over- 
steering car (K < 0) the speed sensitivity of the heading mode is even further exag- 
gerated, while the damping decreases, causing the car to become oscillatory. For an 
undersceer ing car (K > 0) the speed sensitivity of the heading mode denominator is 
reduced and the damping increases with speed. 

Another factor to consider is the heading rate sensitivity to steering inputs. If 
we evaluate the derivative of Equation 6 at zero frequency we end up with simple 
expressions for steady-state heading rate and side acceleration: 
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Hero wo see that steering sensitivity is a function of speed, wheelbase (a 4' b - l), 
and the stability factor, K. At low speeds the car follows a path whose curvature (C) 
IS proportional to wheel deflection and inversely proportional to wneelbase (i.e., 
Ackernann steering): 


C " 6v,/t 


"•60C>- 



Then the heading rate due to following this curved path Is proportional to velocity 


Therefore, 


r - CUq - 

- - 


At higher speeds the stability factor exerts additional speed effects (Equation 13), 
with understeering cars (K > 0) having less speed sensitivity and oversteering cars 
(K < 0) having Increased speed sensitivity. 


Driver Behavior 

Given the above approximate lateral response dynamics let us now analyze the 
corresponding driver control dynamics for the V and Yj. blocks of Figure 1. The Y. 
block Includes several components of driver behi^lor. Farst there Is a component due 
to basic limitations In the driver's response properties. These limitations can be 
subdivided further Into subcomponents: pure time delays due to central nervous system 

processing and neural conduction to the limbs; and interface effects due to the spring 
mass damping system formed by the driver's arms coupled to the steering system. Most 
of the neuromuscular effects are high frequency and can be approximated by a pure time 
delay, e“^®, In the frequency range of Interest for car control (Reference A). 

The second component is a lead or anticipation term (TtS + 1) that the driver 
adopts to counteract vehicle response characteristics discussedi above. The third corn- 

heading 
function 

that has Seen developed and used in a variety of past studies (e.g*. References 2 and 
5): 


ponent is a gain which sets the magnitude of 6 corrections for given 
errors (♦p). Combining the above components we derive a heading response 
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A pure gain 

feedback for lane 

position errors has been found satisfactory 

in pre- 

vious research: 
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The addition of weighted lane position error and heading angle error can actually be 
considered as a composite angular error as shown in Figure 1 : 
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(16) 


It has been previously shown that this equivalent angular error can be Interpreted 
as the angular error to a projected aim point located some distance down the road » 
as shown in Figure 3 (Reference 2). The aim point concept is appealing because it 
represents somewhat of a perceptual efficiency for the driver. Instead of separately 


Figure 


Aim 
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perceiving both lane position and heading errors* he need only perceive the angular 
error to the aim point down the road. 


One additional term will be added here that has received only scant e'^tentlon in 
the Uterature (References 6 and 7). Analysis and modeling of driver response ta 
at ST I has shown the need for low-frequency characteristics or compensation that will 
act to reduce lane position error ofisets and speed up the driver model's transient 
response and error reduction in tasks such as lar. ? changes. The need for low fre- 
quency compensation can be seen by considering the effect of constant input distur- 
bances, 4 j, to the Figure 1 block diagram. Given a constant 6j input the driver model 
(i.e., Yy^and YJ discussed so far would have to allow constant lane position errors 
in order' to generate a compensating wheel angle 6^ at the differential summing block. 
In the real world this would mean that drivers subjected to steady crosswinds or 
crowned roadways (i.e., inputs causing a constant force input to the vehicle) would 
drive with a constant lane position error. 


The above offset error effect is not very 
reasonable, and the model can be corrected to 
eliminate it by adding a parallel trimming inte- 
grator. as illustrated in Figure 4, somewhere in 
the driver model feedforward path. The effect of 
the parallel integrator is to continue increasing 
its output in the face of steady errors, and 
holding its value as the errors approach zero. 
This will then produce steady wheel angles. 6 , 
to compensate for steady disturbances, 6^, wltH- 
out requiring a steady offset error. 


♦u 


^ 


Figure 4. Parallel or "Trim 
Integrator for Counteracting 
Steady- State Error 


There are two possible locations for the parallel Integrator. One is in the \ 
block that would operate only on lane position errors. Another possibility is the 
block, where the parallel Integrator would operate on the composite signal which is 
a combination of lane position and heading errors. The Y^ location seems more reason- 
able for two reasons. First, the concept of perceiving^a simple aim point error 
would still be valid, which would not be true if the parallel Integrator were applied 
to Just the lane position errors. Second, steady-state heading errors can also 
develop in situations such as following curved paths, and the Y. parallel integrator 
location would tend to compensate for these errors quicker than waiting for the head- 
ing error to accumulate into lane position errors which are operated on by the Yy 
block. ^ 


Although we have rationalized the inner-loop (Y. block) as the best location for 
Che parallel integrator characteristic, we will analyze both possible locations (Y. 
and Y^) below in order to accumulate further evidence for the best location. Summar- 
izing* the driver response behavior for both parallel Integrator locations we have: 


Inner-Loop (Y^) Parallel Integrator: 

Ky - Ky 


(16) 


Outer-Loop (Y ) Parallel Integrator. 


Y^ - K^(TlS ♦ 1 ‘e 


Y a 5 .t— K_L, K 

Yy - ^ Ky 




(17) 


- 602 - 




Driver /Vehicle Dynamics 


Given the above dynamic characteristics for the vehicle and driver, let us now 
analyze the overall driver/vehlcle system response given by Equation 3. The first 
term in Equation 5 (Y.) is a driver characteristic, and the third is the car heading 
response (Kg is simply^ the steering ratio). The last term in Equation 5 is a combina- 
tion cf driver and vehicle characteristics. The interaction or these characteristics 
when Ky is a pure gain has been considered previously (Reference 2). Here we will 
consider the case where the parallel integrator is in the outer loop. 

Combining Equations 5, 12, and 17 for the Y parallel integrator we end up with the 
following expression for the bracketed term in Equation 5 


y G| 


bKy(8 f K')[s 2 -I- (g/Tf) (Up/bTr)] 


s^(s + ) 


08 ) 


This can be seen to be a classical root locus problem, with one first-order zero (at 
K'), a second-order zero pair, two poles at the origin and a pole due to the heading 
time constant (s « -1/T ). We can now vary the lane position gain (K^,) and plot the 
Locus of roots for Equation 18. Root locus plots for vehicle characteristics used in 
past research (Reference 8) are compared in Figure 5. Here we see that with the 
parallel integrator in the outer loop a complex pair of low-frequency zeros occurs at 
reasonable values of Ky. 


With Che parallel integrator in the Inner loop (Y 
real zeros, one at K' and the other closed-loop zero 
without the oarallel integrator term. However, as 
that for a large heading time constant Che outer- 
results in a complex zero with fairly low damping, 
low-frequency phase curve has a steep slope. However 
that the low-frequency phase curves are never very s 
two real roots as opposed to the Figure 5 complex roo 


. Equation 16) the result is two 
cue to the Equation 18 expression 
illustrated in Figure 5, we see 
•loop parallel integrator always 
This implies that the driver’s 
, data considered below will show 
tteep and can only be fitted with 
ts . 





Figure 3. Root Locus for Outer-Loop Parallel Integrator 
(Vehicle dynamlci described In Figure 6) 
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Phase Marg,in Approximation 

The vehicle heading response (G|) dictates the high-frequency lead compensation 
CO be provided by the driver. Plots of GJ for various levels of vehicle heading 
time constant (T^^) studied in previous research (Reference 8) are given in Figure 6. 
The numerator zero and second-order denominator combine Co give a net flrst-order- 
appeating transfer function combined with the Icinernatic integration, which gives an 
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Figure 6. Test Vehicle Heading Response Transfer Functions 


idditLonal free s in the denominator. Inspection of Che Figure 6 phase curves shows 
chit "ho heading response dynamics ca.‘ be described by an equivalent time constant, 
-^hich should dictate required driver lead compensation (l.e,. Tl • 

Givtn the above components the composite driver/vehicle open- loop transfer func- 
tion c in be written as follows; 
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Driver/Vehicle Interaction 
Due to Outer-Loop Closure 


where a ■ Uq, This equation can be rearranged according to frequency characteris- 
tics: ^ 
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( 20 ) 


Nov assume Che driver a .justs to compensate for combined phase lag characteristics 
of Terns (3) and @ . The 45 deg phase lag point of @ ^ sonew’nat higher chan T^' 
and Jefine3 , Additional phase lead is derived from @ so that the T£^ frequency 
break may be S(^evhat higher than T^ or , 

Given the complete transfer function expression above , it is now useful to con- 
sider an extended crossover model approximation. Assuir.e that Expressions Q) and 
(?) combine into an equivalent first-order lag which in cancelled by the driver’s lead 
term : 
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Then the dr i .’cr / vehic le equivalent open- loop transfer function reduces to 


where 
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for a neutral steering car and la an effective system time delay that accounts for 

the driver's time delay (t) ana residual phase lags (or lead) left over from the 
approximations in Equation 21. 

Now we can analyze Equation 22 to determine the stability limit for the driver fC. 
gain. For stability the unity magnitude of Equation 22 must occur below the 180 deg 
phase lag point. At high frequencies (« >> K' and a) the unity gain point is given by 
-j (i.e., the unity gain crossover frequency) and the Equation 22 phase is given by 
c 

- can“l — - tan“l ^ (23) 

2 wc wc • 


Phase margin is defined at the gain crossover frequency which occurs at relatively 
high frequency, so that 


tan 



«c 


taP’l — 
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Thus the ohase margin , or amount of extra phase shift allowable before Lnstability Is 
reached, Is given oy 


♦ 
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2f 


5 “ ^e^-c 


(24) 


The second term on the right side of Equation 24 la the phaae lag due to the driver's 
low frequency behavior in controlling lane devlationi. The third term is the phase 
lag duo to the driver's basic time delay limitation, %fhich defines the limiting band- 
width he/she can achieve. Note that the outer- loop operntions (K* and Ky) add phase 
lag and further limit the achievable bandwidth, so that the driver must ^radc off 
inner- and outer- loop gains (i.e., and K*, Ky) In order to optimize perfc ince. 

Model Velidation 

A field study using an instrumented car on a closed test course has been previous- 
ly reported on in Reference 8. The dynamics of the test vehicle could be easily modi- 
fied, and the conditions given in Figure 6 were included in a test program involving 
A male5 and 8 females, ages 25-40, with an average of 13 years driving experience. 
Describing functions were obtained for each set of vehicle dynamics using a 
measureioent technique reported previously (Reference 2). 

Averaged driver describing function data for each set of vehicle dynamics are 
shown in Figure 7, along with curve tits according to the model of Equations 19 and 
20. Some data reinterpretation over that reported in Reference 8 was necessary in 
order to obtain detailed model fits. Although each vehicle data aet was fit indivi- 
dually, Home constraints were observed across vehicles in order to obtain parameter 
values that changed in an orderly manner with vehicle heading time conatanC. 

Model parameters are plotted as a function of inverse equivalent vehicle time con- 
stant (I.e., the basic bandwidth of vehicle heading or yawing response) in Figure 8. 
In general, the trends shown in Figure 8 are consistent with the known behavior of the 
human operator (Reference 4), to wit: 

a Lead generation (Ti) Increases with Increasing system lag (T*q) , 
although the cancellation is not complete. ^ 

a Oporacor ciiae delay (t) increases with lead generation (T|^)< 

a Open- loop gain (A K ) decreases with Increasing system lag to maintain 
stability. ^ 


In this case of mult Iple- loop dynamics we also note that the outer-loop gain ) 
decre.ises with Increasing vehicle lag, which according to Equation 24 also tends Co 
maint.iln stability. 
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Figure 7. Field Teat Drlver/VehicLe Describing Functions and 
Model Fits ( Af‘«malc) 


The equivalent crossover model parametirs plotted in Figure 9 give further insight 
into the driver's adaptation to differe'^t vehicle dynamics. Here ae see that tne 
phase lag component due to the driver’s lane control behavior ((K* + 
naintained relatively constant. This is accomplished by reducing Ky as is reduced 
in response to increased system time delay (t«) which in turn result.) from increased 
vehicle lags. 

CONCLUDltlC REMiUUCS 

The appro<imate dr iver/vehicle steering dynamics analySi^s developed herein pro* 
vides insight into the driver's adaptive behavior. The driver offsets increased 
volucle lags with anticipation or lead behavior, b»»t in doing so incurs additional 
time delay penalty. The driver then compensates for the phase lag due to extra time 
vivlay by reducing his/her gain. 

The effect of the driver's Inherent time delay penalty on stability is analyzed 
with a phase margin approximation. This approximation shows that both crosjover fre- 
pio \c/ and outer-loop gain affect steering stability. Thus, the driver's benavlor in 
controlling lane position results In a phase lag penalty which influences the dlrec- 
ti.’inil stability o: the driver/vehicle system. 

The analysis in this paper relates to directional control stability Independent of 
ciif pith the driver is commanded to follow. Path commands due to roadway curvature 
evi)xe additional driver behavior which has been considered previously in References 3 

ind 3. 
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Figure 8. Effects of Vehicle Dynamics on Driver Model Parameters 



Figure 9. Effects of Vehicle Dynamics on Driver 
Crossover Model Parameters 
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